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The  validation  of  complex  nonlinear  numerical  models  such  as  the  National  Center  for  Atmospheric 
Research  thermospheric  general  circulation  model  (NCAR  TGCM)  requires  a  detailed  comparison  of 
model  predictions  with  data.  The  Equinox  Transition  Study  (ETS)  of  September  17-24.  1984,  provided 
a  unique  opportunity  to  address  the  verification  of  the  NCAR  TGCM,  since  unusually  large  quantities 
of  high-quality  thermospheric  and  ionospheric  data  were  obtained  during  an  intensive  observation 
interval.  In  a  companion  paper  (paper  I)  by  Crowley  et  al.  (this  issue)  a  simulation  of  the  September 
18-19  ETS  interval  was  described.  Using  a  novel  approach  to  modeling,  the  TGCM  inputs  were  tuned 
where  possible  with  guidance  from  data  describing  the  appropriate  input  fields,  and  the  arbitrary 
adjustment  of  input  variables  in  order  to  obtain  thermospheric  predictions  which  match  measurements 
was  avoided,  in  the  present  paper  the  winds,  temperatures,  and  densities  predicted  by  the  TGCM  are 
compared  with  measurements  from  the  ETS  interval.  In  many  respects,  agreement  between  the 
predictions  and  observations  is  good.  The  quiet  day  observations  contain  a  strong  semidiurnal  wind 
variation  which  is  mainly  due  to  upward  propagating  tides.  The  storm  day  wind  behavior  is 
significantly  different  and  includes  a  surge  of  equatorward  winds  due  to  a  global  propagating 
disturbance  associated  with  the  storm  onset.  The  density  data  confirm  the  existence  of  the  newly 
discovered  four-cell  high-latitude  density  anomaly  described  in  paper  I.  A  quantitative  statistical 
comparison  of  the  predicted  and  measured  winds  indicates  that  the  equatorward  winds  in  the  model 
are  weaker  than  the  observed  winds,  particularly  during  storm  times.  This  is  consistent  with  predicted 
latitudinal  temperature  gradients  and  storm  lime  density  increases  which  are  much  smaller  than  the 
observed  values.  Soft  particle  precipitation  or  high-altitude  plasma  healing  is  invoked  as  a  possible 
source  of  the  additional  high-latitude  heating  required  by  the  model.  A  quiet  day  phase  anomaly  in  the 
measured  F  region  winds  which  is  not  reproduced  by  the  model  suggests  the  occurrence  of  an 
important  unmodeled  interaction  between  upward  propagating  semidiurnal  tides  and  high-latitude 
effects.  Wind  data  from  altitudes  below  100  km  indicate  shortcomings  in  the  generic  equinox  solar 
minimum  tidal  specification  used  in  the  TGCM.  The  lack  of  appropriate  data  to  specify  input  fields 
seriously  impairs  our  ability  to  generate  realistic  global  thermospheric  simulations.  The  problem  is 
particularly  acute  in  the  southern  hemisphere.  Furthermore,  disturbances  generated  in  the  southern 
hemisphere  simulation  propagate  northward  and  degrade  the  northern  hemisphere  predictions.  Thus, 
if  the  thermosphere  in  the  northern  hemisphere  is  ever  to  be  fully  understood,  the  southern  hemisphere 
needs  to  be  observed  and  simulated  more  accurately.  Several  improvements  are  suggested  lor  future 
realistic  lime-dependent  simulations  of  specific  intervals. 
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I.  Introduction 

The  validity  of  complex  numerical  models,  such  as  the 
National  Center  for  Atmospheric  Research  (NCAR)  thermo¬ 
spheric  general  circulation  model  (TGCM),  must  be  assessed 
so  that  they  can  be  exploited  for  further  geophysical  studies 
and  for  practical  applications,  such  as  “space  weather" 
prediction.  Validation  of  thermospheric  models  is  particu¬ 
larly  difficult  since  thermospheric  science  is  data-poor,  es¬ 
pecially  when  compared  with  tropospheric  studies.  This  is 
partly  due  to  the  difficulty  of  measuring  basic  parameters 
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such  as  density  and  temperature  at  thermospheric  heights. 
The  Equinox  Transition  Study  (ETS)  of  September  17-24, 
1984,  provided  a  unique  opportunity  to  address  the  verifica¬ 
tion  of  the  NCAR  TGCM,  since  unusually  large  quantities  of 
high-quality  thermospheric  and  ionospheric  data  were  ob¬ 
tained  during  an  intensive  observation  interval  | Carlson  and 
Crowley ,  this  issue). 

One  of  the  main  problems  for  obtaining  realistic  predic¬ 
tions  of  global  thermospheric  dynamics  using  TGCMs  is  the 
specification  of  the  three-dimensional,  time-dependent  forc¬ 
ing  fields  for  any  given  day.  There  are  enough  adjustable 
parameters  in  the  specification  of  TGCM  forcing  fields  that 
the  simulation  of  thermospheric  effects  at  a  single  station  is 
relatively  easy.  However,  little  effort  has  so  far  been  made  to 
reduce  the  number  of  free  parameters  (by  reference  to  large 
quantities  of  data)  and  then  to  accurately  simulate  thermo¬ 
spheric  conditions  at  many  locations  simultaneously. 

In  a  companion  paper  by  Crowley  et  al.  Ithis  issue], 
hereinaftercalled  paper  I,  the  thermospheric  response  to  the 
ETS  magnetic  storm  of  September  19,  1984,  was  simulated 
using  the  NCAR  TGCM.  The  inputs  for  this  simulation  were 
“state  of  the  art,”  and  the  number  of  free  parameters  was 
reduced  to  a  minimum.  The  arbitrary  adjustment  of  input 
variables  in  order  to  obtain  thermospheric  predictions  which 
match  measurements  was  explicitly  avoided. 

The  location  and  size  of  the  auroral  oval  and  the  energies 
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of  precipitating  particles  in  the  model  were  tuned  by  refer¬ 
ence  to  satellite  data  collected  during  the  ETS  period.  The 
resulting  auroral  electron  density  enhancements  were  added 
to  the  densities  obtained  from  the  international  reference 
ionosphere  (IR1)  [Bilitza,  1986],  and  the  total  electron  den¬ 
sity  was  used  to  calculate  ion  drag,  momentum  forcing,  and 
Joule  heating  terms  at  each  TGCM  time  step.  The  0-0 + 
collision  cross  section  is  another  important  parameter  in 
determining  the  ion  drag.  Recent  work  by  Burnside  el  al. 
[1987]  suggests  that  the  0-0 +  collision  cross  section  of 
Schunk  and  Walker  [1973]  is  too  small,  and  its  value  was 
therefore  increased  by  a  factor  of  1.5  in  the  present  simula¬ 
tion.  The  convection  patterns  [Knipp  et  al.,  this  issue]  for 
the  ETS  interval  were  derived  from  the  assimilative  mapping 
of  ionospheric  electrodynamics  (AMIE)  technique  of  Rich¬ 
mond  and  Kamide  [1988],  which  provides  a  best  fit  potential 
pattern  (in  the  least  squares  sense)  to  a  variety  of  iono¬ 
spheric  data.  Paper  l  also  emphasized  the  importance  of 
upward  propagating  semidiurnal  tides  for  thermospheric 
dynamics.  These  tides  were  specified  at  the  TGCM  lower 
boundary  using  the  results  of  Fesen  el  al.  [1986]  for  equinox 
solar  minimum  conditions. 

In  the  present  paper  the  accuracy  of  the  TGCM  predic¬ 
tions  is  assessed  by  comparison  with  observations  from 
many  locations.  Meridional  neutral  winds  are  the  main 
geophysical  parameter  used  in  assessing  the  validity  of  the 
TGCM  simulations  for  a  number  of  reasons:  they  may  be 
inferred  by  a  variety  of  measurement  techniques,  they  are 
available  at  all  times  of  day  from  many  locations,  they  are 
very  sensitive  to  changes  of  high  latitude  and  tidal  forcing, 
and  they,  in  turn,  play  a  major  part  in  determining  temper¬ 
ature,  density,  and  composition  changes  on  a  global  scale. 
For  the  ETS  interval,  meridional  neutral  wind  estimates 
were  obtained  from  the  European  Incoherent  Scatter  (E1S- 
CAT),  Sondre  Stromfjord,  Millstone  Hill,  and  Arecibo  inco¬ 
herent  scatter  radars  in  addition  to  a  network  of  12  ion- 
osondes  in  the  northern  hemisphere  (Figure  1).  For  the  first 
time,  the  differences  between  the  predicted  and  measured 
winds  are  quantitatively  assessed.  Times  and  locations 
where  the  differences  are  likely  to  be  exceptionally  small  or 
large  are  identified. 

In  addition  to  the  meridional  winds,  the  predicted  neutral 
temperatures  and  densities  are  also  compared  with  data.  The 
model  accurately  reproduces  many  features  of  the  ETS  data 
and  may  therefore  be  considered  realistic.  Analysis  of  the 
model  thus  permits  various  observed  phenomena  to  be 
investigated  and  leads  to  an  improved  understanding  of  the 
thermosphere.  Discrepancies  between  the  predictions  and 
observations  yield  additional  insights.  Areas  requiring  fur¬ 
ther  research  are  identified,  and  several  improvements  are 
suggested  for  future  simulations. 

2.  Results 

2. 1 .  Meridional  Neutral  Winds 

The  neutral  wind  velocity  component  in  the  magnetic 
plane  may  be  determined  from  incoherent  scatter  radar 
measurements  of  F  region  ion  drift  [e.g,,  Vusseur,  1969; 
Harper,  1973;  Salah  and  Holt,  1974;  Oliver  and  Salali, 
1988].  The  technique  has  been  widely  applied  at  various 
radar  locations.  It  involves  the  determination  of  the  ion 
velocity  parallel  to  the  magnetic  field  at  a  particular  altitude 


Fig.  1.  Location  of  12  ionosondes  (solid  circles)  and  four  inco¬ 
herent  scatter  radars  (open  circles)  from  which  meridional  neutral 
winds  were  derived  during  the  ETS. 


in  the  F  region,  and  subtraction  of  the  diffusion  velocity 
component.  Ignoring  the  normally  small  vertical  wind  com¬ 
ponent,  the  projection  of  the  residual  velocity  onto  the 
horizontal  plane  yields  an  estimate  of  the  neutral  wind. 

Computation  of  the  ion  diffusion  velocity  involves  a 
knowledge  of  the  neutral  atmospheric  temperature  and  com¬ 
position,  in  addition  to  a  value  for  the  0-0  ’  collision  cross 
section  [e.g.,  Burnside  et  al.,  1988].  The  neutral  parameters 
are  normally  obtained  from  a  model  such  as  MS1S  ( Hedin , 
1983,  1987],  which  may  be  significantly  different  from  the 
conditions  which  actually  existed  on  the  date  of  interest.  The 
0-0 1  collision  cross  section  is  currently  a  subject  of  debate. 
Burnside  et  al.  [1987]  have  suggested  that  the  values  derived 
from  the  Schunk  and  Walker  [1973]  formula  should  be 
increased  by  a  factor  of  1.7  [  +®  j], 
lonosonde  data  have  also  been  used  to  deduce  the  merid¬ 
ional  neutral  wind  component.  Miller  et  al.  \  1986]  describe  a 
method  for  obtaining  the  wind  from  measurements  of  the 
height  of  the  maximum  electron  concentration  in  the  Fz 
region  (/imjx).  The  method  involves  the  use  of  an  ionospheric 
model  to  develop  a  relation-Eir  'setween  the  meridional  wind 
speed  and  /imax.  The  mer i  !-<..nat  vind  is  then  derived  by  a 
comparison  of  the  modele.’  measured  layer  heights. 
Miller  et  al.  [1987]  estimate  i..  it  statistical  errors  in  the 
ionosonde  data  analysis  generate  wind  uncertainties  of  30%. 
This  value  was  used  to  represent  enor  bars  in  the  present 
paper,  although  more  realistic  error  bars  might  be  as  large  as 
75  m  s' 1  when  some  of  the  assumptions  involved  in  the 
analysis  are  considered.  The  effects  of  vertical  winds  on  hmdX 
are  neglected,  as  in  the  radar  analysis.  The  derivation  of 
neutral  winds  by  this  technique  also  neglects  the  effects  of 
electric  fields,  although  Miller  et  al.  [1987]  found  that  this 
does  not  usually  lead  to  significant  errors  during  quiet 
conditions.  The  analysis  also  involves  many  of  the 

same  assumptions  as  the  radar  analysis  in  terms  of  neutral 
atmospheric  parameters  and  the  0-0  *  collision  cross  sec¬ 
tion.  The  assumptions  and  uncertainties  involved  in  obtain¬ 
ing  meridional  neutral  winds  from  both  incoherent  scatter 
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radar  ion  drift  data  and  data  are  discussed  in  detail  by 
Buonsanto  et  al.  [this  issue]. 

Figure  2  illustrates  the  meridional  neutral  winds  (geomag¬ 
netic  north  positive)  as  a  function  of  universal  time  measured 
(solid  line)  at  five  locations  during  September  18,  a  quiet  day, 
and  September  19,  a  disturbed  day.  The  TGCM  predictions 
(dashed  line)  are  also  plotted  for  comparison.  Millstone  Hill 
and  Arecibo  are  incoherent  scatter  radar  sites  in  the  Amer¬ 
ican  sector,  while  Rome,  Tomsk,  and  Tokyo  are  ionosonde 
sites  in  the  European,  Russian,  and  Japanese  sectors,  re¬ 
spectively  (Table  1).  Rome  and  Millstone  Hill  have  similar 
geographic  latitudes,  but  the  magnetic  latitude  of  Rome  is 
about  15°  less  than  that  at  Millstone  Hill.  Arecibo  and  Tokyo 
have  the  same  magnetic  latitude,  but  the  geographic  latitude 
of  Arecibo  is  about  20°  less  than  that  of  Tokyo.  Note  that  the 
uncertainties  in  the  winds  derived  from  ionosonde  data  are 
generally  comparable  with  those  associated  with  incoherent 
scatter  radar  data,  and  the  predictions  for  each  station  fall 
within  the  measurement  error  bars  at  various  times,  espe¬ 
cially  on  the  quiet  day. 

The  quiet  day  measurements  display  a  strong  diurnal 
variation,  with  a  tendency  toward  poleward  daytime  winds. 
Table  1  indicates  the  universal  and  local  times  for  which 
quiet  time  poleward  and  equatorward  winds  were  observed 
for  each  station.  The  duration  of  the  poleward  winds  is 
shorter  for  higher  magnetic  latitudes.  Thus  Millstone  Hill 
observed  poleward  winds  for  only  9  hours,  between  05  and 
13  LT  on  the  quiet  day,  in  contrast  to  the  other  stations 
where  winds  were  poleward  for  up  to  18  hours.  The  equa¬ 
torward  winds  of  the  quiet  day  reached  magnitudes  of 
100-150  m  s~'  at  each  station. 

The  large  poleward  winds  obtained  from  the  Rome  and 
Tokyo  ionosonde  data  in  the  evening  hours  (Figures  2  c  and 
2e)  were  unexpected,  and  their  origin  appears  to  be  tidal. 
The  data  do  not  appear  to  be  erroneous;  however,  the 
analysis  neglected  electric  field  effects  which  may  have 
contributed  to  the  apparent  poleward  winds. 

The  observed  strong  diurnal  wind  variation  is  reasonably 
well  modeled  by  the  TGCM.  An  additional  feature  of  the 
predicted  winds  is  the  strong  semidiurnal  variation  which 
was  discussed  in  paper  I.  Maxima  in  the  model  poleward 
winds  at  about  06  and  18  LT  are  generated,  in  part,  by  the 
upward  propagating  semidiurnal  tides.  As  a  result  of  these 
tides,  the  model  predicts  poleward  winds  during  about  16 
hours  for  each  site.  This  agrees  well  with  the  observed 
variation  listed  in  Table  I  for  Arecibo,  Tokyo,  Rome,  and 
Tomsk,  although  the  model  wind  variation  tends  to  be  phase 
shifted  to  earlier  times  (0-4  hours)  with  respect  to  the  data. 
For  Millstone  Hill  the  predicted  winds  contain  a  strong 
semidiurnal  component  which  contrasts  with  the  observed 
wind  variation.  The  predicted  winds  thus  blow  poleward  for 
16  hours,  compared  with  the  observed  9  hours  which  is 
typical  of  summertime  behavior  for  this  location. 

On  the  active  day  the  measured  winds  are  considerably 
different  from  their  quiet  day  values.  At  Millstone  Hill  the 
winds  remained  equatorward  during  the  whole  day,  diverg¬ 
ing  from  the  quiet  day  wind  pattern  as  early  as  05  UT.  This 
behavior  suggests  that  some  form  of  disturbance  occurred  to 
the  north  of  Millstone  Hill  before  the  main  storm  com¬ 
menced  at  10  UT.  Sondre  Stromfjord  (not  shown)  is  the  only 
other  site  containing  similar  early  divergence  from  quiet  day 
values,  indicating  that  the  “precursor”  event  was  localized 
in  the  American  sector. 


Fig.  2.  Meridional  neutral  winds  close  to  300  km  measured 
(solid)  and  predicted  (dashed)  for  (a)  Millstone  Hill,  ( b )  Arecibo,  (c) 
Rome,  ( d )  Tomsk,  and  ( e )  Tokyo.  Negative  values  indicate  equa¬ 
torward  winds;  units  are  meters  per  second.  Note  that  the  radar  data 
(Millvone  Hill  and  Arecibo)  were  averaged  over  50  km  in  altitude 
and  smoothed  with  a  l-hour  running  mean  before  interpolating  to 
hourly  values.  Vertical  arrows  indicate  local  midnight. 
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TABLE  !.  Intervals  During  Which  Measured  Winds  Plotted  in  Figure  2  Are  Poleward  or  Equatorward 


Geographic 

Coordinates 

Magnetic  Coordinates 

Declination 

Quiet  Time  Poleward 
Winds 

Quiet  Time  Equatorward 
Winds 

Latitude 

Longitude 

Latitude 

Longitude 

UT 

LT 

UT 

LT 

Millstone 

42.6 

-71.5 

53.9 

-1.4 

-14.0 

10-18 

05-13 

19-09 

14-04 

Arecibo 

18.3 

-66.7 

29.8 

3.9 

-8.0 

09-24 

05-20 

01-08 

21-04 

Tokyo 

36.4 

140.6 

28.5 

-150.7 

-7.0 

18-11 

03-20 

12-17 

21-02 

Rome 

41.9 

12.52 

36.05 

86.5 

-3.0 

06-20 

07-21 

21-05 

22-06 

Tomsk 

56.5 

84.9 

46.0 

160.6 

12.0 

00-13 

04-17 

14-23 

18-03 

The  most  outstanding  feature  of  the  winds  measured  on 
September  19  is  the  surge  of  equatorward  winds  between 
about  10  and  17  UT.  This  surge  was  present  in  every  data  set 
and  was  associated  with  a  large  disturbance  which  propa¬ 
gated  equatorward  from  the  auroral  zone.  The  surge  will  be 
discussed  more  fully  in  a  later  paper,  but  we  briefly  mention 
some  of  the  surge  features  here,  in  comparing  the  measured 
and  predicted  winds. 

Winds  derived  from  ionosonde  data  were  available  only 
hourly.  The  radar  data  were  therefore  averaged  and  interpo¬ 
lated  to  hourly  values  and  compared  with  hourly  TGCM 
predictions.  The  magnitude  of  the  observed  wind  surge  in 
Figure  2  is  100-200  ms-1,  with  peak  velocities  occurring 
later  at  lower  latitudes.  Analysis  of  the  unaveraged  data  (not 
shown)  reveals  that  the  disturbance  has  two  maxima,  both  of 
which  reached  Arecibo  1  hour  after  passing  the  location  of 
Millstone  Hill.  This  corresponds  to  a  propagation  speed  of 
650-700  m  s'.  At  the  present  time,  the  interpretation  of  the 
two  maxima  in  the  data  is  unclear.  The  two  peaks  may  be 
consecutive  cycles  of  the  same  wave,  or  they  may  be  due  to 
entirely  independent  waves.  For  the  purposes  of  the  present 
paper  the  distinction  is  relatively  unimportant,  since  the 
TGCM  reproduces  the  first  maximum  very  well,  at  all 
locations,  but  the  second  wind  maximum  is  not  predicted. 
The  predicted  surge  therefore  appears  to  last  for  a  shorter 
time  than  the  observed  event. 

Figure  3  is  a  geographic  polar  plot  depicting  the  progres¬ 
sion  of  the  peak  of  the  disturbance  from  high  to  low  latitudes 
between  12  and  14  UT  on  September  19.  The  TGCM 
predictions  are  denoted  by  solid  lines,  while  the  maxima 
deduced  from  the  wind  data  are  given  by  dashed  lines.  The 
predicted  phase  fronts  are  distorted  because  of  different 
propagation  speeds  on  the  dayside  and  nightside.  Part  of  the 
distortion  probably  arises  from  propagation  through  the 
strong  antisunward  polar  cap  wind  system.  The  source  of  the 
disturbance  in  the  TGCM  is  the  postmidnight  sector,  and  this 
is  the  sector  in  which  the  wave  amplitudes  are  largest  and  the 
phase  progression  is  clearest.  The  amplitude  of  the  distur¬ 
bance  is  much  smaller  on  the  opposite  side  of  the  polar  cap 
in  the  afternoon  sector,  and  the  shape  of  the  phase  fronts 
cannot  be  traced  with  certainty.  Figure  2  c  indicates  that  the 
peak  of  the  predicted  surge  at  Rome  occurred  about  13  UT 
and  had  an  amplitude  of  only  75  m  s-1.  The  situation  is 
further  confused  in  this  sector  because  the  phase  fronts 
appear  to  be  discontinuous.  At  90°E  the  disturbance  is 
present  simultaneously  at  two  different  latitudes,  presum¬ 
ably  because  the  energy  took  different  paths. 

Because  of  the  incomplete  data  set  and  I -hour  time 
resolution,  the  sketching  of  measured  wave  fronts  in  Figure 
3  is  incomplete  and  somewhat  ambiguous;  however,  the 


model  appears  to  have  predicted  many  features  of  the 
observed  disturbance  with  remarkable  accuracy.  Globally, 
the  distorted  wave  fronts  are  reproduced,  with  a  discontinu¬ 
ity  in  the  European  sector.  The  speed  of  propagation  is 
virtually  identical  with  the  measured  speed,  and  the  timing  of 
the  event  is  correct  to  within  15  min. 

The  storm  time  wind  variation  in  the  Japanese  and  Russian 
sectors  is  exceptionally  complex,  as  illustrated  by  Figures  2e 
and  2d,  for  Tokyo  and  Tomsk,  respectively.  At  Tokyo  the 
propagating  disturbance  is  superposed  on  the  normal  mid¬ 
night  surge,  such  that  the  latter  appears  to  be  increased  in 
amplitude.  However,  the  model  winds  quickly  abate  because 
the  second  part  of  the  surge  was  not  generated  in  the  TGCM 
and  also  because  of  the  arrival  of  a  disturbance  generated  in 
the  southern  hemisphere.  We  surmise  that  the  high-latitude 
forcing  in  this  sector  of  the  southern  hemisphere  may  be  too 
strong  in  the  model. 

The  winds  observed  at  Tomsk  were  enhanced  throughout 
the  afternoon  and  evening  from  12  to  23  UT  beginning  with 
the  storm-associated  surge  which  occurred  about  12  UT  (18 
LT).  In  the  model  the  surge  is  predicted  at  12-13  UT;  then 
the  winds  return  to  quiet  time  values  before  further  strong 
equatorward  winds  are  predicted  from  16  to  21  UT.  The 


O* 


Fig.  3.  Global  propagation  of  disturbance  between  12  and  14  UT 
on  September  19  (day  263).  Latitude  contours  are  0°,  30°,  and  60°N. 
Solid  curves  represent  the  locus  of  maximum  equatorward  winds 
and  temperature  increases  predicted  by  the  TGCM,  Dashed  lines 
illustrate  wind  maxima  deduced  from  data.  Solid  circles  indicate 
locations  from  which  wind  data  were  available. 
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discrepancy  between  the  observed  and  predicted  winds  from 
14  to  16  UT  is  again  due  to  both  the  lack  of  a  second  “cycle” 
in  the  predicted  storm  surge  and  the  disturbance  which  in  the 
model  propagates  from  the  southern  hemisphere.  A  similar 
picture  is  obtained  throughout  this  longitude  sector. 

It  is  interesting  to  note  that  the  generation  of  large 
disturbances  in  the  southern  hemisphere  of  the  TGCM 
occurred  near  local  midnight  for  both  the  Japanese  and 
European  sectors.  A  second  disturbance  launched  at  18  UT 
in  the  southern  hemisphere  reached  Rome  at  22-23  UT, 
causing  a  rapid  abatement  of  the  model  winds,  while  the  data 
indicate  an  unusually  strong  equatorward  wind  at  that  time 
(Figure  2c).  The  southern  hemisphere  high-latitude  forcing 
in  the  midnight  local  time  sector  may  be  generally  too  strong 
in  the  TGCM. 

Exact  replication  of  measured  thermospheric  properties 
by  the  TGCM  should  not  be  expected,  since  both  the  data 
analysis  techniques  and  the  model  are  imperfect.  The  veri¬ 
fication  of  thermospheric  models  and  the  intercomparison  of 
various  simulations  require  a  quantitative  analysis  of  the 
differences  between  the  measured  and  predicted  properties, 
or  between  model  predictions.  We  therefore  present  a  brief 
statistical  analysis  of  the  differences  between  the  measured 
winds  of  the  ETS  period  and  those  predicted  by  the  TGCM 
simulation.  This  analysis  has  been  performed  for  the  com¬ 
bined  48-hour  data  set,  applying  equal  weighting  to  data  from 
the  four  incoherent  scatter  radars  and  the  12  ionosondes 
(Figure  1).  In  general,  the  rms  and  mean  differences  will  be 
discussed,  instead  of  the  standard  deviation.  This  is  because 
we  are  interested  in  how  well  the  TGCM  could  predict 
thermospheric  properties  in  a  forecasting  situation.  Compu¬ 
tation  of  the  standard  deviation  requires  a  subtraction  of  the 
mean  difference,  which  would  not  be  known  a  priori. 

Figure  4  indicates  the  distribution  of  the  differences  be¬ 
tween  measured  and  predicted  winds  for  day  262  and  day  263 
with  their  respective  mean  and  rms  deviations.  The  differ¬ 
ences  range  from  +250  ms"1  to  -350  m  s"1,  but  both  the 
mean  and  rms  deviations  are  smaller  on  day  262  (September 
18)  than  on  day  263  (September  19).  For  the  active  day  a 
mean  deviation  of  -27  ms-1  indicates  that  the  observed 
equatorward  winds  are  consistently  stronger  than  those 
produced  by  the  model. 

The  distribution  of  differences  has  been  broken  down 
further  into  smaller  sets,  in  order  to  identify  under  what 
circumstances  they  are  most  significant.  The  rms  difference 
has  different  properties  on  days  262  and  263  (Figure  5).  On 
the  quiet  day  (262)  it  is  approximately  constant  with  latitude, 
while  for  the  active  day  (263)  the  rms  difference  increases 
with  latitude.  These  variations  are  particularly  clear  when 
plotted  in  terms  of  geographic  latitude.  The  mean  difference 
(not  shown)  becomes  more  negative  with  increasing  latitude 
on  both  days  and  is  on  average  30  ms'1  more  negative  on 
day  263  at  all  latitudes.  The  same  result  is  obtained  for  both 
geographic  and  geomagnetic  latitudes. 

The  rms,  mean,  and  standard  deviations  of  the  differences 
are  plotted  as  functions  of  universal  time  in  Figure  6.  The 
number  of  stations  available  for  each  UT  is  also  presented 
(Figure  6 d)  and  varies  from  10  to  15.  The  rms  deviation 
increases  from  about  60  m  s  " 1  on  the  quiet  day  to  more  than 
100  m  s"1  after  the  storm  on  September  19.  The  rms 
deviation  comprises  the  mean  and  the  standard  deviations, 
both  of  which  also  increase  over  48  hours.  The  most  dra¬ 
matic  increase  is  in  the  mean  deviation,  which  is  approxi¬ 


A  WIND  (m*-1) 

Fig.  4.  Distribution  of  differences  between  measured  and  pre¬ 
dicted  meridional  neutral  wind  speeds.  The  mean  (0.2  ms"1)  and 
rms  (59  m  s'1)  deviations  for  the  quiet  day  (262)  are  smaller  than 
those  for  the  active  day,  which  were  -27  m  s"1  and  87  m  s~l. 
respectively. 


mately  zero  before  the  storm  onset  and  during  the  equator- 
ward  surge  ( 10—13  UT,  day  263)  but  then  rapidly  drops  to 
almost  - 100  m  s  ' 1  before  slowly  recovering  to  -50  ms'1. 
These  negative  values  indicate  that  the  measured  winds 


GGLAT 

Fig.  5.  Latitudinal  variation  for  rms  deviation  of  the  predicted 
winds  from  the  measured  winds,  (u)  Quiet  day  (262).  ( b )  Active  day 

(263). 
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Fig.  6.  Statistical  deviations  of  predicted  winds  from  measured  winds  for  48  hours  during  September  18-19,  1984:  (a) 
rms  difference,  (b)  mean  difference,  (c)  standard  deviation,  and  (d)  number  of  stations  contributing  to  statistics. 


remained  strongly  equatorwards  after  the  storm  for  longer 
than  the  TGCM  predicts. 

Dividing  the  stations  into  three  approximately  equal 
groups  at  high,  middle,  and  low  latitudes,  and  considering 
each  day  separately,  the  local  time  variation  of  the  differ¬ 
ences  is  examined  in  Figure  7.  For  the  quiet  day  (Figure  la) 
the  mean  difference  generally  lies  between  ±50  m  s  at  all 
latitudes,  with  a  tendency  to  be  largest  for  middle  and  high 


latitudes  at  night.  For  low  latitudes  the  mean  difference 
tends  to  maximize  around  dawn  and  dusk  and  may  be 
associated  with  the  phasing  of  the  tides  in  the  TGCM,  as 
discussed  for  individual  stations  above. 

On  the  active  day  (Figure  lb)  the  differences  are  smallest 
at  low  latitudes  and  largest  at  high  latitudes.  The  same 
distinct  diurnal  variation  of  the  mean  differences  occurs  at  all 
three  latitude  ranges;  the  smallest  differences  occur  during 


Fig.  7.  Local  lime  variation  of  mean  difference  between  predicted  and  measured  winds.  Low,  middle,  and  high 
geographic  latitudes  were  considered  separately,  (u)  Quiet  day  (262).  (b)  Active  day  (263).  Negative  values  indicate  that 
the  measured  winds  blew  more  strongly  equatorward  than  the  predicted  winds. 
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Fig.  8.  Exospheric  temperature  variation  predicted  and  measured 
for  Arecibo  during  September  18-19,  1984. 


the  daytime  between  08  and  18  LT.  A  similar  conclusion  is 
reached  for  the  rms  and  standard  deviations  (not  shown). 

2.2.  Neutral  Temperatures 

The  determination  of  exospheric  and  neutral  temperatures 
from  incoherent  scatter  radar  measurements  has  recently 
been  described  in  detail  for  Arecibo  [ Burnside  et  al.,  1988] 
and  for  Millstone  Hill  [Hunan  and  Salah,  1988;  Oliver  and 
Salah,  1988],  In  this  section  the  exospheric  temperature 
estimates  from  the  Arecibo  and  Millstone  Hill  ETS  data  are 
compared  with  those  predicted  by  the  TGCM  and  MSIS-86 
models. 

Figure  8  illustrates  the  exospheric  temperature  variations 
predicted  for  September  18-19,  1984,  at  Arecibo  by  the 
MSIS-86  model  and  the  TGCM.  Both  models  predict  mod¬ 
erate  temperature  increases  of  about  50  K.  over  quiet  time 
values  following  the  storm  onset  at  10  UT  on  September  19. 
In  general,  the  TGCM  predictions  exceed  the  MSIS  temper¬ 
atures  by  about  50  K.  Incoherent  scatter  radar  data  were 
available  from  10  UT  on  September  19  and  are  superposed 
on  Figure  8.  The  measured  temperatures  are  about  50  K  less 
than  those  predicted  by  the  MSIS  model,  and  almost  100  K 
less  than  those  of  the  TGCM  except  for  the  last  few  hours  of 
September  19.  The  data  also  contain  two  short-lived  temper¬ 
ature  enhancements  at  about  13  UT  and  !5UTon  September 
19,  which  are  associated  with  the  equatorward  wind  surge 
described  above.  As  with  the  winds,  the  TGCM  predicted 
the  first  peak  (13  UT),  but  not  the  second  (15  UT). 

A  similar  discrepancy  between  the  MSIS  model  and 
measured  temperatures  was  reported  both  by  Hunan  and 
Salah  [1988]  and  by  Burnside  et  al.  [  1988]  for  solar  minimum 
conditions.  Burnside  et  al.  (1988]  attribute  the  discrepancy 
to  a  problem  with  the  MSIS  model.  Analysis  of  the  Arecibo 
radar  data  incorporated  in  the  MSIS  model  neglected  H  * 
ions  for  altitudes  around  400  km,  since  they  were  not 
originally  thought  to  be  significant  at  solar  minimum.  This 
gave  rise  to  MSIS  temperatures  systematically  too  hot  by  50 
K.  During  the  ETS  interval,  solar  F  10.7  fluxes  were  compa¬ 
rable  with  those  of  solar  minimum  conditions  (about  70),  and 
significant  concentrations  of  H +  were  observed  at  Arecibo 
[ Tepley  and  Kerr,  1987].  Temperatures  deduced  from  the 
Arecibo  radar  data  were  not  corrected  for  H +  and  will 


Fig.  9.  Exospheric  temperature  variation  predicted  and  measured 
for  Millstone  Hill  during  September  18-19,  1984. 


therefore  be  too  high  by  about  50  K.  Note  that  the  discrep¬ 
ancy  between  the  MS1S/TGCM  predictions  and  measured 
temperature  would  have  been  even  larger  if  H  *  ions  had 
been  taken  into  account  in  the  ETS  data  analysis. 

Figure  9  illustrates  the  corresponding  temperature  varia¬ 
tions  for  Millstone  Hill.  The  TGCM  predictions  essentially 
track  MSIS,  except  after  21  UT  on  September  19  when  the 
TGCM  recovers  more  rapidly  from  storm  effects  than  MSIS. 
For  both  nights  the  quiet  time  mode)  temperatures  are  100  K. 
hotter  than  the  measured  temperatures.  Conversely,  the 
models  underestimate  the  observed  storm  effects  by  up  to 
350  K,  subject  to  the  uncertainties  discussed  below.  A 
similar  discrepancy  was  noted  by  Hunan  [1988]  for  a  winter 
solstice  solar  minimum  study. 

The  validity  of  storm  time  exospheric  temperature  esti¬ 
mates  from  radar  observations  has  been  discussed  in  detail 
by  Hanan  [1988],  who  suggests  that  overestimates  may 
result  from  contamination  of  the  ion  temperatures  by  ion- 
frictional  heating,  particle  bombardment,  and  molecular 
ions.  In  the  present  case,  none  of  these  effects  is  strong 
enough  to  generate  the  exospheric  temperature  enhance¬ 
ments  plotted  in  Figure  9.  Both  the  frictional  heating  and 
molecular  ion  contamination  would  require  electric  fields  of 
about  50  mV  m " 1  and  must  be  discounted  because  the 
electric  fields  measured  at  Millstone  Hill  were  too  small  (<20 
mV  m  "').  Similarly,  although  there  was  evidence  of  particle 
precipitation  in  the  vicinity  of  Millstone  Hill  during  the  ETS 
interval,  it  was  not  sufficiently  intense  to  generate  the 
observed  exospheric  temperature  increase.  These  increases 
will  be  discussed  in  detail  in  a  separate  paper  and  are  due  to 
two  causes.  The  background  temperature  increases  of  100- 
200  K  between  12  and  24  UT  are  due  to  the  transport  of  hot 
air  from  the  auroral  zone  into  the  Millstone  Hill  radar  field  of 
view.  An  additional  enhancement  at  14  UT  is  caused  by 
compression  associated  with  the  second  peak  of  the  propa¬ 
gating  wind  surge  which  was  not  predicted  by  the  TGCM. 

The  temperature  gradient  between  Millstone  Hill  and 
Arecibo  is  of  interest  because  it  affects  the  meridional 
neutral  wind  component.  For  the  quiet  day  the  two  models 
predict  temperature  differences  of  ±75  K,  depending  on 
local  time.  Both  models  predict  that  on  the  active  day  the 
Millstone  Hill  temperatures  increase  more  than  the  Arecibo 
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temperatures,  resulting  in  a  temperature  gradient  which 
tends  to  drive  stronger  equatorward  winds.  The  MSIS  model 
predicts  that  the  Millstone  Hill  temperature  remains  50  K 
higher  than  the  Arecibo  temperature,  while  the  TGCM 
predicts  a  smaller  difference.  In  contrast,  the  measured 
temperature  difference  between  Arecibo  and  Millstone  Hill 
is  typically  200  K,  reaches  a  maximum  of  450  K,  and  falls  to 
75  K  by  24  UT.  In  addition,  the  Arecibo  temperatures  are 
probably  overestimated  by  ~50  K.  owing  to  the  neglect  of 
H  +  ,  as  mentioned  above,  so  the  true  temperature  gradient 
may  be  even  larger.  The  radar  data  from  September  19  are 
thus  consistent  with  the  models  in  the  sense  that  the  Mill¬ 
stone  Hill  temperatures  increased  more  than  the  Arecibo 
temperatures;  however,  the  magnitude  of  the  observed  in¬ 
crease  and  corresponding  temperature  gradient  is  much 
greater  than  the  models  predict.  The  temperature  data  are 
thus  qualitatively  consistent  with  stronger  equatorward 
winds  than  predicted  by  the  models,  especially  during  the 
storm  period. 

Rishbeth  [1972]  provides  formulae  to  estimate  the  neutral 
winds  that  would  result  from  a  given  horizontal  F  region 
temperature  gradient.  A  simplistic  calculation  indicates  that 
the  75-K  to  400-K  temperature  difference  between  Millstone 
Hill  and  Arecibo  would  drive  meridional  neutral  winds 
ranging  from  150  ms-1  to  820  ms1.  It  is  therefore  possible 
to  account  quantitatively  for  the  lack  of  equatorward  winds 
in  the  TGCM  as  compared  with  the  data. 

2.3.  Neutral  Density 

During  the  ETS  interval  the  Air  Force  Geophysics  Labo¬ 
ratory  (AFGL)  S85-I  satellite  accelerometer  experiment 
[ Marcos  and  Forbes,  1985]  measured  total  atmospheric 
density  in  an  almost  circular  orbit  at  altitudes  between  197 
and  201  km.  Examination  of  the  TGCM  predictions  for  this 
altitude  region  has  led  to  the  discovery  of  a  new  high-latitude 
density  anomaly  as  described  briefly  in  paper  I.  The  mor¬ 
phology  of  the  anomaly  and  its  causative  mechanisms  will  be 
discussed  in  detail  in  a  later  paper. 

Several  aspects  of  the  satellite  data  merit  inclusion  in  the 
TGCM  verification  process.  Figure  10«  illustrates  the  mea¬ 
sured  density  variation  with  latitude  in  the  noon  sector  at 
about  1420  UT.  The  typical  quiet  day  variation  is  given  by  a 
dashed  line  for  comparison.  By  1420  UT  the  density  at  200 
km  has  increased  by  30-50%  at  all  latitudes  because  of  the 
storm  time  energy  input  and  resultant  thermospheric  expan¬ 
sion.  Two  high-density  regions  were  observed  at  65°N  and 
40°N. 

To  set  these  observations  in  context,  the  TGCM  density 
predictions  for  200  km  at  14  UT  are  presented  as  Figure  10b. 
In  this  case,  the  dominant  polar  cap  features  are  a  high- 
density  region  in  the  noon  sector  and  an  intense  low-density 
region  in  the  dawn  sector.  The  satellite  path  is  traced  as  a 
heavy  line  in  the  noon  sector  and  indicates  that  the  high 
density  observed  at  65°N  is  part  of  the  high-latitude  density 
anomaly  noted  above,  while  the  density  peak  centered  at 
40°N  is  due  to  an  equatorward  propagating  disturbance. 
When  the  altitude  variation  of  the  satellite  is  taken  into 
account,  the  amplitudes  of  the  observed  and  predicted 
density  perturbations  are  also  very  similar. 

Figure  1 1  compares  the  typical  variation  of  observed 
prestorm  densities  with  those  from  the  late  storm  phase 
about  17-21  UT.  The  densities  measured  at  200  km  during 
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Fig.  10.  (a)  Latitudinal  variation  of  neutral  density  measured  by 

the  S85-I  satellite  for  altitudes  close  to  200  km  about  142 0  UT  on 
September  19,  1984.  The  typical  corresponding  quiet  day  variation  is 
also  shown  for  comparison  (dashed  line).  (/>)  Northern  hemisphere 
density  (x It)-15  g  cm'5)  at  200  km  from  the  TGCM  simulation  for 
14  UT  on  September  19,  1984.  The  satellite  track  from  Figure  I0«  is 
superimposed  on  the  dayside. 


the  storm  have  increased  by  50%  on  both  the  nightside  and 
the  dayside.  In  contrast  the  densities  in  the  TGCM  (not 
shown)  increased  by  only  about  16%.  The  relatively  small 
increase  in  the  TGCM  densities  indicates  that  insufficient 
energy  has  been  deposited  into  the  TGCM  during  the  simu¬ 
lation. 

3.  Discussion 

The  final  state  predicted  by  any  model  depends  to  some 
extent  on  the  initial  state.  In  the  present  case,  the  TGCM 
was  run  for  several  days  from  a  diurnally  reproducible 
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Fig.  II.  Latitude  variation  of  neutral  density  measured  close  to 
200  km  by  the  S85-I  satellite.  Data  from  passes  at  1730  UT  (solid 
line)  and  2130  UT  (dashed  line)  on  September  19  are  compared  with 
the  typical  corresponding  quiet  day  variation. 

equinox  solar  minimum  steady  state,  using  inputs  appropri¬ 
ate  for  a  quiet  day  during  September  1984.  Since  magnetic 
conditions  were  quiet  for  several  days  preceding  the  ETS 
interval  and  since  we  have  concentrated  here  on  winds, 
which  respond  fairly  rapidly  to  changes  in  the  model  inputs 
[Mayr  and  Hedin ,  1977],  the  initial  conditions  are  not 
expected  to  cause  serious  problems  for  the  present  study. 
This  assumption  is  borne  out  by  the  relatively  good  agree¬ 
ment  between  the  measured  and  predicted  winds  and  tem¬ 
peratures  for  the  quiet  day  (September  18)  simulation. 

In  order  to  assess  the  validity  of  the  TGCM  simulation  for 
the  ETS  interval,  the  predicted  winds,  temperatures,  and 
densities  have  been  compared  with  those  measured  during 
the  ETS  campaign.  In  many  respects,  the  agreement  be¬ 
tween  the  models  and  measurements  is  excellent,  but  several 
discrepancies  have  also  emerged  from  this  analysis.  The 
main  features  will  now  be  discussed. 

3.1.  Quiet  Day  Winds 

The  quiet  day  F  region  wind  variations  predicted  by  the 
model  display  a  strong  semidiurnal  variation  similar  to  that 
observed  at  many  locations,  indicating  that  tidal  e fleets  are 
extremely  important  and  must  be  modeled  correctly  if  ther¬ 
mospheric  predictions  are  to  be  realistic.  In  the  present  case, 
the  phase  of  the  modeled  quiet  day  wind  variation  '.eads  the 
measured  phase  by  0-4  hours,  and  this  discrepancy  needs  to 
be  addressed.  The  problem  is  particularly  clear  in  the 
nighttime  equatorward  winds,  where  the  largest  phase  shift 
seems  to  be  at  Rome,  while  the  smallest  shift  is  at  the  higher 
magnetic  latitude  stations:  Tomsk  and  Millstone  Hill.  A 
similar  problem  has  been  noted  for  other  intervals  (e.g.,  V. 
B.  Wickwar,  Global  thermospheric  studies  of  neutral  dy¬ 
namics  using  incoherent  scatter  radars,  submitted  to  Ad¬ 
vances  in  Space  Research,  1989),  which  suggests  that  the 
problem  is  significant.  The  local  time  variation  of  the  Mill¬ 
stone  predictions  leads  the  predicted  variation  for  Rome  by 
about  I  hour.  This  would  be  expected  because  of  the 
clockwise  diurnal  wind  rotation  and  the  different  magnetic 
declinations  at  the  two  sites  (Table  I)  if  the  two  wind 
systems  were  identical  in  geographic  cooruinates.  The  data 
reveal  that  the  measured  wind  variation  at  Millstone  Hill 


leads  the  Rome  variation  by  about  4  hours,  which  is  more 
difficult  to  explain. 

The  phase  of  the  thermospheric  winds  is  determined  by 
the  interaction  of  various  diurnal  and  semidiurnal  compo¬ 
nents.  The  thermospheric  diurnal  tides  are  excited  almost 
exclusively  in  situ  [Lindzen,  1970;  Lindzen  and  B'uke,  1971] 
except  for  the  equatorial  lower  thermosphere.  In  this  study, 
only  in  situ  diurnal  excitation  is  included.  Semidiurnal  tidal 
oscillations  are  generated  by  (1)  an  in  situ  component  due  to 
solar  UV  and  EUV  absorption,  (2)  a  propagating  component 
due  to  the  upward  moving  waves  excited  by  insolation 
absorption  in  the  lower  atmosphere,  and  (3)  diurnal  varia¬ 
tions  in  ion  drag  interacting  with  diurnal  variations  in  neutral 
velocities.  The  in  situ  excitation  and  ion-neutral  momentum 
coupling  are  explicitly  included  in  the  model:  the  propagat¬ 
ing  semidiurnal  components  at  the  lower  boundary  are 
parameterized.  Verification  of  the  in  situ  tidal  forcing  for  the 
ETS  interval  is  beyond  the  scope  of  this  paper;  however,  the 
propagating  semidiurnal  component  at  the  lower  boundary  is 
amenable  to  investigation  and  is  discussed  below. 

Since  Millstone  Hill  and  Rome  have  the  same  geographic 
latitude,  the  observed  phase  anomaly  cannot  easily  be  ex¬ 
plained  purely  in  terms  of  upward  propagating  tides,  al¬ 
though  Herrero  et  al.  [1983]  have  qualitatively  explained 
longitudinal  differences  in  tidal  behavior  by  invoking  a 
coupling  between  the  upward  propagating  tides  and  the  tides 
generated  in  situ  by  both  EUV  heating  and  the  interaction  of 
the  diurnal  wind  field  variations  with  the  diurnal  ion  density 
variations.  The  predicted  phase  anomaly  appears  to  be  fairly 
insensitive  to  ion  drag,  and  therefore  the  background  iono¬ 
sphere  in  the  TGCM.  Numerical  experiments  revealed  that 
increasing  the  ion  drag  by  50%  on  a  global  scale  produced 
insignificant  changes  in  the  relative  phases  of  the  tidal  winds 
predicted  for  F  region  altitudes.  Since  Millstone  Hill  has  a 
higher  magnetic  latitude  than  Rome,  differences  in  the  ef¬ 
fects  of  high-latitude  forcing  at  the  two  sites  appear  to  offer 
the  simplest  explanation  of  the  observed  phase  anomaly. 

Recent  work  by  C.  G.  Fesen  (private  communication, 
1988)  indicates  that  the  amplitude  and  phase  of  the  semidi¬ 
urnal  tide  obtained  in  the  F  region  are  very  sensitive  to  the 
degree  of  high-latitude  forcing.  High-latitude  forcing  can 
generate  a  pseudotide  that  interacts  with  the  stable  tide  from 
below,  and  a  Fourier  decomposition  of  the  disturbed  state 
gives  an  apparent  phase  shift  in  the  resultant  semidiurnal 
component.  The  phase  of  this  component  can  change  by  up 
to  3  hours  at  mid-latitudes  as  the  cross-polar  cap  potential 
drop  is  increased  from  zero  to  90  kV.  The  amplitude  of  the 
apparent  semidiurnal  tide  may  correspondingly  increase  by  a 
factor  of  2.  The  discrepancy  between  the  modeled  and 
observed  phases  of  the  F  region  winds  in  the  present  study 
may  well  be  attributable  to  inaccurately  simulated  interac¬ 
tions  of  this  type  in  the  TGCM.  High-latitude  effects  were 
not  considered  in  the  analysis  of  Fesen  et  al.  [1986],  where 
the  equinoctial  tide  specification  was  originally  deduced. 
Clearly,  the  level  of  geomagnetic  activity  should  be  consid¬ 
ered  in  future  studies  when  deriving  generic  specifications 
for  upward  propagating  tides  by  matching  F  region  data  with 
model  simulations.  Furthermore,  the  generic  tidal  specifica¬ 
tions  deduced  from  this  type  of  analysis  need  to  be  used  with 
caution. 

For  the  ETS  simulation  the  propagating  semidiurnal  tidal 
inputs  at  the  lower  boundary  were  those  obtained  by  Fesen 
et  al.  [1986]  for  equinox  solar  minimum  conditions.  These 
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values  were  used  because  the  Fl0  7  fluxes  observed  during 
September  1984  were  comparable  with  those  typical  of  solar 
minimum  conditions.  Fesen  et  al.  [1986]  '  sumed  that  the 
antisymmetric  modes  would  be  unimportant  at  the  equinox. 
They  obtained  amplitudes  and  phases  of  the  2,  2  and  2,  4 
modes  as  follows:  A  (2,  2)  =  600  m,  A  (2,  4)  =  300  m,  <t>  (2, 
2)  =  -0.5  hour,  <f>  (2,  4)  =  -2.0  hours.  Richmond  and  Rob/e 
[1987]  investigated  the  electrodynamic  effects  of  thermo¬ 
spheric  winds  from  the  NCAR  TGCM  and  found  that  the 
inclusion  of  this  tidal  forcing  in  the  equinoctial  TGCM 
improved  the  agreement  between  calculated  and  observed 
magnetic  variations  and  electric  fields. 

Since  the  higher-order  modes  tend  to  be  attenuated  below 
300  km,  F  region  winds  are  dominated  by  the  interaction  of 
the  in  situ  tides  with  the  upward  propagating  (2,  2)  mode.  By 
way  of  a  numerical  experiment  we  discovered  that  when  the 
phase  of  the  (2,  2)  mode  is  changed  at  the  TGCM  lower 
boundary,  the  phase  of  the  predicted  F  region  winds  changes 
by  a  corresponding  amount.  The  incorrect  specification  of 
any  of  the  higher  tidal  modes  at  the  TGCM  lower  boundary 
should  not  necessarily  degrade  the  wind  predictions  at  F 
region  heights,  although  they  clearly  influence  the  quality  of 
the  predictions  at  lower  altitudes. 

For  September  1984,  wind  data  were  available  from  alti¬ 
tudes  around  97  km  at  four  locations.  It  was  therefore 
possible  to  compare  the  semidiurnal  tidal  amplitudes  and 
phases  used  at  the  TGCM  lower  boundary  with  the  mea¬ 
sured  monthly  mean  values.  The  data  for  Saskatoon  were 
published  by  Munson  et  al.  [1986]  for  Poker  Flat  by  Teten- 
baum  et  al.  [1986],  and  for  Kyoto  by  Vincent  et  al.  [1988]. 
Adelaide  data  were  generously  provided  by  R.  A.  Vincent 
(private  communication,  1988).  Figure  12  summarizes  the 
measured  phases  and  amplitudes  of  the  eastward  and  north¬ 
ward  semidiurnal  components.  The  phase  variation  with 
altitude  (oblique  solid  line)  was  approximately  the  same  for 
three  stations  (Saskatoon,  Poker  Flat,  and  Adelaide).  For 
Kyoto,  data  were  only  published  for  90  km,  and  the  altitude 
variation  is  assumed  to  be  an  average  of  that  observed  at  the 
other  stations.  The  measured  amplitudes  are  only  indicated 
at  one  height  for  each  station  since  the  variation  with  altitude 
was  generally  small. 

The  phases  of  the  meridional  components  are  essentially 
the  same  for  all  four  sites;  however,  the  northern  hemisphere 
zonal  components  are  approximately  in  antiphase  with  the 
Adelaide  data.  Vincent  et  al.  1 1988]  found  a  similar  relation¬ 
ship  to  be  present  throughout  the  year  when  they  compared 
semidiurnal  tides  from  the  geographically  conjugate  loca¬ 
tions  of  Kyoto  and  Adelaide.  They  also  discovered  that  the 
phases  exhibit  a  bimodal  character  and  change  rapidly  from 
one  state  to  the  other  at  the  equinoxes.  The  zonal  phase  for 
90  km  at  Kyoto  changes  around  September/October  from  a 
value  of  0700  LT  in  summer  to  a  value  near  0000  LT  in 
winter.  The  data  in  Figure  12  thus  indicate  that  the  semidi¬ 
urnal  tides  during  the  ETS  campaign  were  in  a  “northern 
hemisphere  summer"  state.  Vincent  et  al.  [I988|  argue  that 
this  type  of  north-south  asymmetry  can  only  be  explained  by 
the  dominance  of  antisymmetric  modes  in  the  semidiurnal 
wind  oscillations.  They  invoke  the  presence  of  the  (2,  3) 
mode  in  particular,  since  the  (2,  5)  mode  has  nulls  in  the 
horizontal  wind  structure  function  near  35°N. 

Superimposed  on  Figure  12  are  crosses  corresponding  to 
the  phases  and  amplitudes  of  the  simulated  winds  at  the 
TGCM  lower  boundary.  These  winds  correspond  to  the 


upward  propagating  semidiurnal  tides  specified  according  to 
Fesen  et  al.  [1986],  For  the  northern  hemisphere  the  mea¬ 
sured  phase  leads  the  predicted  phase  by  about  1 .0-1 .5  hours 
in  both  components.  This  represents  exceptionally  good 
agreement.  However,  in  the  southern  hemisphere  the  mea- 
suied  and  predicted  winds  are  approximately  in  antiphase 
for  both  components.  The  winds  at  the  lower  boundary  of 
the  TGCM  thus  fail  to  reproduce  the  observed  hemispheric 
asymmetry,  and  this  is  probably  because  the  antisymmetric 
modes  were  excluded  from  the  TGCM  tidal  input  specifica¬ 
tion.  The  phase  problem  in  the  F region  is  more  complicated, 
however.  We  note  that  the  F  region  and  97-km  phase 
differences  between  the  TGCM  and  observed  winds  in  the 
northern  hemisphere  have  the  opposite  sign.  Thus  a  simple 
correction  to  the  total  phase  at  97  km  would  exacerbate  the 
F  region  problem. 

In  addition  to  the  phase  problem,  the  amplitudes  of  the 
different  modes  need  to  be  specified.  For  each  site  in  Figure 
12  the  amplitude  of  the  TGCM  winds  at  97  km  exceeds  the 
measured  amplitude  by  factors  ranging  from  30 %  to  200%. 
Richmond  and  Ruble  [1987]  have  also  noted  that  the  semid¬ 
iurnal  tidal  amplitudes  as  specified  by  Fesen  et  al.  [1986] 
should  be  reduced  by  about  30%.  A  reduction  in  the  total 
semidiurnal  amplitude  at  97  km  may  arise  from  introducing 
antisymmetric  modes  if  their  phase  is  appropriate  for  de¬ 
structive  interference.  A  detailed  analysis  of  the  ETS  wind 
data  from  altitudes  below  100  km  in  order  to  improve  the 
tidal  specification  at  the  TGCM  lower  boundary  is  outside 
the  scope  of  this  paper.  However,  this  type  of  analysis  may 
prove  very  fruitful  in  the  future,  particularly  if  data  could  be 
obtained  from  more  sites.  The  Lower  Thermosphere  Cou¬ 
pling  Study  (LTCS),  a  receri  CEDAR  campaign,  may  be 
able  to  address  this  issue. 

3.2.  Active  Day  Winds 

For  the  storm  day  the  “state  of  the  art"  TGCM  simulation 
has  predicted  many  features  of  the  observed  winds  very 
well.  Particularly  noteworthy  is  the  reproduction  of  the 
discontinuous  wave  fronts  displayed  in  Figure  3.  following 
the  storm  onset.  The  disturbance  also  has  the  correct  ampli¬ 
tude  and  propagation  speed,  which  suggests  that  the  high- 
latitude  inputs  are  very  realistic  at  this  time.  However, 
comparison  of  the  predictions  with  data  has  revealed  several 
discrepancies  in  the  model  winds:  (I)  there  is  no  “precur¬ 
sor"  event  at  06  UT  in  the  American  sector,  (2)  the  second 
part  of  the  obser/ed  equatorward  wind  surge  (-14  UT)  is 
missing,  (3)  disturbances  generated  in  the  southern  hemi¬ 
sphere  have  amplitudes  which  are  too  large  when  they  reach 
the  northern  hemisphere,  and  (4)  northern  hemisphere  equa¬ 
torward  winds  are  too  weak  later  in  the  storm  period, 
particularly  during  the  nighttime  at  high  latitudes. 

The  precursor  event  of  September  19  noted  in  connection 
with  Figure  2a  is  a  departure  from  the  quiet  day  wind 
variation  at  Millstone  Hill  at  06  UT,  prior  to  the  onset  of  the 
magnetic  storm.  This  disturbance  appears  to  be  confined  to 
the  high-latitude  American  sector,  including  Sondre  Strom- 
fjord  and  Millstone  Hill.  The  AF  index  plotted  in  paper  I 
indicates  that  a  small  substorm  commenced  about  05  UT  on 
September  19  and  peaked  at  06  UT  with  a  value  of  300  nT. 
Magnetograms  for  24  stations  throughout  the  northern  hemi¬ 
sphere  were  available  during  this  interval  | Knipp  et  al.,  this 
issue].  They  revealed  that  activity  was  confined  to  the 
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Fig.  12.  Monthly  mean  (September  1984)  phase  and  amplitude  variation  as  a  function  of  altitude  for  eastward  and 
northward  wind  components  at  (a)  Poker  Flat,  (6)  Saskatoon,  (c)  Kyoto,  and  (J)  Adelaide.  Oblique  lines  represent 
measured  phase  variation  with  altitude.  Solid  circles  represent  measured  values;  crosses  indicate  values  in  TGCM. 


midnight  sector,  corresponding  to  longitudes  between 
Greenland  and  Churchill  (50°-100°W).  Unfortunately,  no 
magnetograms  were  available  for  the  immense  region  be¬ 
tween  Greenland  and  Churchill  and  north  of  Ottawa  (55°N 
magnetic),  although  the  Millstone  Hill  radar  measured  elec¬ 
tric  fields  of  less  than  10  mV  m"1  up  to  65°N  (magnetic). 

The  convection  patterns  used  in  the  TGCM  were  derived 
from  the  AMIE  technique  and  depend  on  good  data  coverage 
for  their  accuracy.  The  variations  of  cross-polar  cap  poten¬ 
tial  and  Joule  heating  obtained  from  the  AMIE  procedure 
were  described  in  paper  I  and  by  Knipp  et  al.  [this  issue). 
Both  parameters  have  small  values  at  06  UT.  The  potentials 
and  E  fields  obtained  from  AMIE  during  this  time  were 
underestimated,  since  the  model  had  essentially  no  informa¬ 
tion  about  the  magnetic  activity  in  the  critical  sector.  It  is 
interesting  to  note  that  the  winds  measured  at  Millstone  Hill 


after  06  UT  on  September  19  continued  to  blow  equatorward 
throughout  the  prestorm  period.  The  reason  for  this  is  still 
not  clear,  since  substorm  activity  seems  to  have  ceased 
between  0630  and  0930  UT. 

The  equatorward  wind  “surge"  observed  globally  be¬ 
tween  10  and  17  UT  contained  two  peaks  separated  by  about 
2  hours  and  will  be  discussed  in  more  detail  in  a  later  paper. 
The  nature  of  this  surge  is  not  clear,  and  the  two  current 
working  hypotheses  regard  the  surge  either  as  two  cycles  of 
a  single  periodic  wave  or  as  two  distinct  events  each  having 
a  single  cycle.  The  surge  generated  in  the  TGCM  simulation 
corresponded  to  the  first  of  the  observed  peaks.  The  source 
of  the  disturbance  in  the  TGCM  was  a  region  of  intense  Joule 
healing  over  Alaska  and  western  Canada  in  the  postmidnight 
sector  at  about  10  UT.  This  was  approximately  the  time  at 
which  the  AE  index  increased  dramatically  by  1000  n'P 
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because  of  the  storm  onset.  The  origin  of  a  second  observed 
disturbance  might  be  expected  to  be  in  a  similar  location  to 
the  source  of  the  first  event  (since  many  of  the  properties  of 
the  two  disturbances  are  similar)  except  2  hours  later  (i.e., 
about  12  UT);  however,  there  is  no  evidence  of  such  a 
source  in  the  high-latitude  data.  The  AE  index  is  generally 
decreasing  from  11  UT.  At  12  UT  there  is  a  short-lived 
enhancement  from  800  to  10  nT,  but  it  seems  unlikely  that 
such  an  event  could  launch  a  disturbance  of  the  same 
magnitude  as  the  10  UT  storm  onset.  This  evidence  indicates 
that  the  second  peak  may  be  the  second  cycle  of  a  single 
large-scale  gravity  wave  having  a  period  of  about  2  hours. 
The  second  cycle  of  the  F  region  wind  variation  was  not 
predicted  by  the  TGCM,  and  theory  predicts  that  dissipation 
should  restrict  the  observed  mid-latitude  response  to  a  single 
pulse  [Testud,  1973;  Richmond,  1978).  However,  the  TGCM 
density  predictions  from  lower  altitudes  did  contain  a  prop¬ 
agating  disturbance  at  14  UT  (Figure  10),  and  it  is  not  clear 
at  the  moment  either  why  there  was  no  corresponding  F 
region  wind  effect  or  what  was  the  source  of  the  disturbance. 

Wavelike  disturbances  have  rarely  been  observed  on  the 
global  scale.  A  recent  example  from  the  Worldwide  Atmo¬ 
spheric  Gravity  Wave  Study  (WAGS)  campaign  was  de¬ 
scribed  by  Rice  et  al.  [1988],  who  deduced  the  most  likely 
origin  of  their  wave  to  be  in  the  evening  sector  of  the  auroral 
oval.  Identification  of  the  source  characteristics  is  a  primary 
goal  of  gravity  wave  research.  Further  analysis  of  the  TGCM 
for  the  ETS  period  promises  to  reveal  important  information 
regarding  the  source  of  the  propagating  disturbance(s)  de¬ 
scribed  here. 

lonosonde  data  from  the  southern  hemisphere  indicate 
that  a  large-amplitude  disturbance  was  also  generated  at 
southern  high  latitudes  in  association  with  the  storm  onset 
about  10  UT.  The  TGCM  predicted  a  large-scale  disturbance 
launched  in  the  southern  hemisphere  about  10  UT  which 
propagated  equatorward.  The  amplitude  of  the  model  distur¬ 
bance  was  largest  in  the  Australian  sector  (~I50°E),  al¬ 
though  it  had  only  half  the  amplitude  of  the  corresponding 
disturbance  in  the  northern  hemisphere.  Both  disturbances 
crossed  the  equator  at  about  the  same  time  (14  UT)  and 
propagated  into  the  opposite  hemisphere.  The  disturbance 
from  the  southern  hemisphere  reached  Tokyo  at  16  UT, 
causing  a  rapid  reversal  from  equatorward  to  poleward 
winds  (Figure  2e)  in  the  TGCM  predictions.  A  similar 
problem  was  encountered  in  the  European  sector  about  22 
UT  (Figure  2c)  due  to  a  southern  hemisphere  disturbance 
launched  at  18  UT  in  association  with  the  third  major 
substorm  onset  of  the  day.  These  disturbances  from  the 
southern  hemisphere  thus  contribute  to  the  discrepancy 
between  the  predicted  and  measured  winds  in  the  northern 
hemisphere. 

As  explained  in  paper  l,  the  lack  of  southern  hemisphere 
high-latitude  data  precludes  the  direct  determination  of 
potential  patterns  using  the  AMIE  technique.  Furthermore, 
the  AMIE  potentials  obtained  from  the  northern  hemisphere 
cannot  be  used  directly  in  the  southern  hemisphere.  As  an 
interim  measure,  the  southern  hemisphere  convection  pat¬ 
tern  was  approximated  by  the  Heelis  model  appropriately 
tuned  by  means  of  the  northern  hemisphere  AMIE  poten¬ 
tials.  The  integrated  Joule  heating  and  momentum  forcing  in 
the  southern  hemisphere  is  thus  expected  to  be  of  the  correct 
order  of  magnitude;  however,  the  TGCM  is  not  expected  to 
reproduce  details  of  the  energy  and  momentum  inputs. 


Therefore  further  analysis  of  the  waves  emanating  from  the 
southern  hemisphere  is  not  deemed  worthwhile  at  this  time. 
However,  the  existence  of  these  waves  indicates  that  the 
southern  hemisphere  cannot  be  ignored  in  the  future  if 
realistic  global  simulations  are  to  be  achieved. 

In  addition  to  the  equatorward  wind  surge  in  the  early  part 
of  the  storm  interval,  unusually  strong  equatorward  winds 
were  also  observed  in  the  late  storm  period  16-23  UT  at 
many  mid-  to  high-latitude  stations.  In  general,  these  winds 
were  underestimated  by  the  TGCM,  particularly  for  the 
nightside.  The  techniques  by  which  meridional  winds  are 
estimated  from  radar  and  ionosonde  data  make  several 
assumptions  which  are  likely  to  be  inappropriate  during 
magnetically  active  intervals,  particularly  at  high  latitudes.  It 
seems  unlikely,  however,  that  the  assumptions  could  ac¬ 
count  entirely  for  the  discrepancy,  which  often  amounts  to 
more  than  100  ms-1  for  individual  stations  (Figure  6b).  The 
TGCM  was  therefore  examined  for  causes  of  underestimated 
equatorward  winds. 

Meridional  winds  in  the  nighttime  sector  can  be  generated 
by  both  momentum  and  heat  sources.  Killeen  and  Roble 
[1984,  1986]  demonstrated  that  the  F  region  winds  are 
controlled  by  the  balance  between  ion  drag  and  pressure 
forces,  and  the  midnight  surge  is  caused  by  acceleration  of 
air  parcels  over  the  polar  cap  by  momentum  forcing.  During 
the  entire  ETS  storm  interval  of  September  19,  antisunward 
winds  in  the  polar  cap  exceeded  500  ms1  in  the  TGCM, 
spilling  out  into  the  nightside  sector  at  lower  latitudes  and 
with  reduced  speeds.  Various  high-latitude  heat  sources  also 
contributed  to  these  winds  by  means  of  pressure  gradients 
directed  from  high  to  low  latitudes. 

Between  18  and  22  UT  the  E1SCAT  radar  was  in  the  late 
evening  sector  and  measured  the  antisunward  winds  flowing 
out  of  the  polar  cap.  The  observed  meridional  wind  speed 
ranged  from  300  to  400  ms'1  while  the  TGCM  predicted 
speeds  of  only  150-200  ms'1.  This  comparison  suggests  that 
the  momentum  forcing  in  the  polar  cap  is  not  accurately 
specified  in  the  TGCM.  Both  the  electric  fields  and  the 
electron  concentration  contribute  directly  to  the  momentum 
forcing,  and  both  involve  a  large  degree  of  uncertainty.  The 
polar  cap  potentials  from  the  AMIE  technique  are  uncertain 
because  of  the  lack  of  polar  cap  electrodynamic  data,  and  in 
future  work  of  this  nature,  polar  cap  drifts  such  as  those 
measured  by  digital  ionosondes  [e.g.,  Reinisch  et  al.,  1987] 
should  be  included  in  the  AMIE  analysis.  The  electron 
concentration  in  the  TGCM  is  specified  by  the  1R1  model 
together  with  cusp,  polar  rain,  and  auroral  precipitation.  In 
the  future,  the  NCAR  thermosphere/ionosphere  general  cir¬ 
culation  model  (TIGCM)  of  Roble  et  al.  [I988«]  will  be 
employed,  which  includes  an  interactive  ionosphere.  One  of 
the  features  of  the  interactive  ionosphere  is  the  generation  of 
enhanced  polar  cap  electron  concentrations  |e.g.,  Fuller- 
Rowell  et  al.,  1987],  which  leads  to  the  possibility  of  stronger 
forcing  of  the  neutral  winds  as  they  are  accelerated  across 
the  polar  cap. 

It  is  difficult  to  explain  the  weak  equatorward  winds 
entirely  in  terms  of  momentum  forcing.  The  weak  latitudinal 
temperature  gradients  and  small  density  increases  in  the 
TGCM  also  suggest  that  additional  high-latitude  heating  is 
required.  The  additional  heating  may  be  supplied  by  Joule 
heating,  particle  heating,  plasma  healing  from  wave-particle 
interactions,  or  enhanced  magnetospheric  heat  flows.  The 
Joule  heating  already  in  the  model  seems  likely  to  be 
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reasonably  accurate  since  the  auroral  convection  and  E 
region  precipitation  have  been  tuned  by  reference  to  mea¬ 
sured  quantities.  Smith  et  til.  [1982]  showed  that  moderate 
fluxes  of  soft  ( —200  eV)  particles  at  high  latitudes  could 
generate  significant  equatorward  winds  at  mid-latitudes. 
Since  the  only  form  of  soft  particles  in  the  TGCM  is 
currently  the  cusp  and  drizzle  precipitation  [Roble  and 
Ridley,  1987],  there  is  good  reason  to  place  additional  soft 
particle  precipitation  in  the  auroral  zone.  In  addition  to 
stronger  equatorward  winds,  such  a  source  would  also 
generate  increased  latitudinal  temperature  gradients  in  the  F 
region,  and  particles  deposited  at  150-180  km  would  produce 
additional  neutral  density  enhancements  at  200  km.  Unde¬ 
tected  soft  particle  precipitation  to  the  north  of  Millstone 
Hill  may  also  explain  the  “precursor  event"  at  06  UT  noted 
in  Figure  2  a. 

4.  Summary  and  Conclusions 

This  is  the  third  in  a  series  of  time-dependent  simulations 
for  specific  intervals  [cf.  Ruble  et  al.,  1987,  1988b],  In 
previous  studies,  the  high-latitude  TGCM  inputs  have  been 
parameterized  in  a  semiempirical  fashion,  which  has  led  to  a 
type  of  quasi-averaged  result.  The  TGCM  inputs  in  the 
present  case  represent  the  current  "state  of  the  art.”  The 
high-latitude  inputs  have  been  tuned  so  that  they  more 
closely  resemble  measured  quantities  such  as  auroral  parti¬ 
cle  precipitation  and  high-latitude  electric  fields.  The  data 
have  enabled  the  convection  pattern  and  auroral  oval  to  be 
defined  with  unprecedented  precision  for  this  simulation; 
however,  the  dynamic  and  nonuniform  nature  of  the  auroral 
oval  and  convection  pattern  means  that  they  cannot  be 
defined  exactly  with  the  present  limited  data  set.  Tuning  of 
the  high-latitude  inputs  was  very  labor  intensive,  although 
the  process  has  now  been  developed  to  a  semiautomatic 
level.  Upward  propagating  semidiurnal  tides  have  also  been 
included  in  the  model,  according  to  the  specifications  of 
Fesen  et  at.  [1986]. 

The  model  predictions  for  September  18-19,  1984,  have 
been  compared  with  global  observations  of  the  meridional 
neutral  wind  in  addition  to  neutral  temperature  and  density 
measurements  on  a  more  limited  scale.  In  many  respects,  the 
agreement  between  the  predictions  and  observations  is  gen¬ 
erally  good,  although  various  discrepancies  remain.  High¬ 
lights  include  the  following: 

1 .  The  strong  semidiurnal  quiet  day  F  region  wind  vari¬ 
ation  is  reproduced  by  the  model  and  emphasizes  the  impor¬ 
tance  of  tidal  effects  for  the  dynamics  of  the  upper  thermo¬ 
sphere. 

2.  The  storm  day  equatorward  wind  surge  with  its  dis¬ 
torted  and  discontinuous  fronts  is  reproduced  with  the 
correct  amplitude  and  propagation  speed.  The  timing  of  the 
event  is  correct  to  within  15  min. 

3.  A  four-cell  high-latitude  neutral  density  anomaly  has 
been  discovered.  Corresponding  features  were  measured  by 
satellite  at  200  km  altitude. 

4.  Enhanced  temperatures  were  measured  at  high  lati¬ 
tudes  following  the  storm  onset.  Similar  enhancements  were 
predicted  by  the  TGCM,  permitting  the  cause  of  the  temper¬ 
ature  increase  to  be  investigated. 

For  the  first  time  the  differences  between  the  predicted 
and  measured  winds  have  been  quantitatively  assessed  on  a 
statistical  basis.  The  differences  tend  to  be  largest  during  the 


night  at  high  latitudes  following  a  magnetic  storm.  The 
smallest  differences  occur  during  the  daytime  at  low  lati¬ 
tudes,  particularly  for  magnetically  quiet  conditions.  In  the 
present  study,  the  best  agreement  between  the  predicted  and 
measured  winds  was  found  for  Arecibo. 

Exact  replication  of  measured  thermospheric  conditions 
by  the  TGCM  should  not  be  expected,  even  assuming  all  the 
important  processes  are  included  in  the  model,  since  the 
initial  conditions  and  many  details  of  the  thermospheric 
forcing  are  not  accurately  known.  The  wind  measurements 
should  also  be  treated  with  caution  since  they  involve 
several  assumptions  which  are  likely  to  be  inappropriate 
during  magnetically  active  conditions,  particularly  at  high 
latitudes. 

Several  discrepancies  between  the  predicted  and  observed 
winds  were  discussed  in  detail: 

1.  The  quiet  day  semidiurnal  variation  of  the  F  region 
winds  is  dominated  by  the  upward  propagating  (2,  2)  mode. 
An  anomalous  phase  relationship  between  the  winds  ob¬ 
served  at  different  locations  is  not  reproduced  by  the  model, 
suggesting  that  interactions  between  the  tides  and  the  high- 
latitude  forcing  are  important  and  require  further  study. 

2.  On  the  active  day  a  precursor  event  observed  in  the 
American  sector  was  not  predicted,  because  of  the  lack  of 
input  data. 

3.  The  nature  of  the  equatorward  wind  surge  observed  at 
Millstone  Hill  and  Arecibo  following  the  storm  onset  is  still 
not  fully  understood.  Detailed  analysis  of  the  wind  variation 
has  revealed  two  maxima  separated  by  2  hours;  however, 
there  is  no  strong  evidence  for  a  source  capable  of  generating 
a  second  disturbance  2  hours  after  the  first  event.  This 
suggests  that  the  surge  may  be  a  single  wave  with  two 
cycles.  The  second  peak  was  not  predicted  by  the  model, 
and  this  discrepancy  is  not  understood. 

4.  Throughout  the  storm  period  the  measured  equator- 
ward  winds  were  stronger  than  the  predicted  winds.  One 
reason  for  this  may  be  insufficient  momentum  forcing  at  high 
latitudes  due  to  inappropriate  polar  cap  convection  and 
electron  concentrations.  A  second  reason  was  the  influence 
of  poleward  winds  associated  with  waves  propagating  from 
the  southern  hemisphere.  Since  very  few  data  were  available 
to  define  the  high-latitude  inputs,  the  southern  hemisphere 
was  not  discussed  in  detail.  This  study  indicates  that  the 
southern  hemisphere  cannot  be  ignored  if  the  northern 
hemisphere  thermosphere  is  ever  to  be  fully  understood. 

5.  Both  the  observed  temperature  gradients  and  storm 
time  density  enhancements  indicate  that  insufficient  high- 
latitude  heating  is  being  applied  to  the  model.  This  is 
quantitatively  consistent  with  the  lack  of  equatorward  winds 
in  the  model.  We  have  suggested  that  the  TGCM  requires  an 
additional  source  of  high-latitude  heat  in  the  form  of  soft 
(<200  eV)  particle  precipitation,  but  the  influence  of  other 
high-altitude  plasma  heal  sources  cannot  be  discounted. 

6.  Wind  data  from  altitudes  below  100  km  revealed 
shortcomings  of  the  generic  equinox  solar  minimum  tidal 
specification  used  in  the  TGCM.  Contrary  to  the  usual 
assumptions,  the  asymmetric  modes  are  still  important  at 
equinox.  In  the  future,  wind  data  from  altitudes  below  100 
km  should  be  used  to  improve  the  TGCM  tidal  inputs  for 
specific  campaign  intervals. 

Eventually,  the  TGCM  will  be  used  as  an  operational 
predictive  model.  The  present  study  reveals  the  type  of  input 
information  which  will  be  necessary  to  make  realistic  fore- 
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casts  of  thermospheric  behavior.  For  magnetically  quiet 
intervals  at  low  latitudes  the  solar  and  tidal  forcing  terms  are 
dominant,  and  the  definition  of  high-latitude  TGCM  inputs  is 
relatively  unimportant.  The  tidal  inputs  were  discussed  in 
detail  above,  and  several  problems  were  identified.  The  solar 
EUV/UV  fluxes  are  also  largely  unknown.  They  are  param¬ 
eterized  in  the  TGCM  by  means  of  the  daily  Fw  7  flux  but  in 
practice  are  highly  variable  and  not  currently  measured. 

To  obtain  accurate  thermospheric  forecasts  for  high  lati¬ 
tudes  and  for  geomagnetically  active  conditions,  a  knowl¬ 
edge  of  the  high-latitude  forcing  terms  is  vital.  Of  particular 
importance  are  the  convection  pattern  and  auroral  particle 
precipitation.  The  results  discussed  above  indicate  that  the 
accurate  specification  of  large-scale  features  of  the  auroral 
oval  such  as  location,  size,  and  particle  energies  is  feasible 
using  a  simple  model  with  small  quantities  of  data.  At  this 
stage,  the  relative  importance  of  smaller-scale  features  such 
as  discrete  auroral  arcs  is  not  known.  The  high-latitude 
convection  pattern  seems  to  be  a  less  tractable  problem. 
This  difficulty  with  the  convection  pattern  is  presumably  a 
reflection  of  the  complexity  of  the  convection  as  compared 
with  the  auroral  oval.  Simultaneous  measurement  of  the 
entire  convection  pattern  has  to  be  a  high  priority  in  the 
future.  The  conductances  in  the  TGCM  are  different  from 
those  used  by  the  AMIE  technique  and  are  not  therefore 
electrically  consistent  with  the  potential  pattern  from  AMIE. 
Future  objectives  should  include  electrodynamic  coupling 
between  the  potential  patterns  and  conductances  in  the 
TGCM. 

Development  of  the  TGCM  is  continuing.  Robte  et  al. 
[  1 988a  1  have  recently  described  a  new  version  of  the  model 
(TIGCM)  which  includes  a  fully  interactive  ionosphere  and 
additional  high-latitude  heating  due  to  soft  particle  fluxes. 
The  aeronomical  scheme  in  the  TIGCM  is  also  revised  and 
includes  radiative  cooling  at  5  ^m  due  to  NO,  in  addition  to 
the  l5-/xm  C02  cooling.  The  lack  of  an  interactive  iono¬ 
sphere  has  been  considered  as  a  possible  contributor  to  the 
differences  between  model  predictions  and  measurements  in 
the  present  study.  With  an  interactive  ionosphere  the  con¬ 
vection  of  dayside  F  region  plasma  over  the  polar  cap  would 
produce  higher  polar  cap  electron  concentrations  than  the 
1RI  model.  The  forcing  of  cross-polar  cap  winds  may  there¬ 
fore  be  enhanced  in  an  interactive  model.  The  interactive 
ionosphere  would  not  significantly  affect  the  tidal  wind 
phases. 

Owing  to  the  known  shortcomings  of  the  current  TGCM, 
the  composition  and  densities  have  not  been  discussed  in 
depth  in  the  present  paper.  Moreover,  the  only  composition 
estimates  available  for  day  262  are  those  from  the  Millstone 
Hill  radar.  These  will  be  discussed  at  a  later  stage,  when  the 
new  TIGCM  is  fully  operational. 

This  study  demonstrates  the  value  of  coordinated  global 
campaigns,  especially  those  which  last  for  several  days. 
Direct  measurements  of  the  thermospheric  neutral  winds 
from  Fabry-Perot  interferometers  (FP1)  would  generally 
complement  the  data  from  other  instruments;  however,  no 
suitable  FPI  data  were  obtained  during  the  ETS  campaign. 
Estimates  of  the  meridional  neutral  wind  obtained  from 
h ma*fr2  measured  by  ionosondes  are  also  particularly  useful 
in  extending  the  network  of  wind  observations.  Although  a 
large  number  (>50)  of  ionosondes  are  currently  operational 
worldwide,  estimates  of  h^Fi,  or  the  parameters  (e.g., 
M(3000)F2)  from  which  it  may  be  calculated,  are  not  rou¬ 


tinely  available  from  most  stations.  If  hourly  values  of 
/imaxF2  could  be  obtained  during  the  coordinated  campaigns, 
the  value  of  ionospheric  soundings  would  be  greatly  en¬ 
hanced.  and  the  data  would  play  an  important  role  in 
improving  our  understanding  of  the  thermosphere. 
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